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Abstract 

We calculate the multiple scattering effect on single hadron production in semi-inclusive lepton- 
nucleus deeply inelastic scattering. We show that the quantum interference of multiple scattering 
amplitudes leads to suppression in hadron productions. At the leading power in medium length, the 
suppression can be approximately expressed in terms of a shift in z of the fragmentation function 
D(z), and could be therefore interpreted as the collisional energy loss. We compare our calculation 
with existing experimental data. We also discuss the effect of quark mass on the suppression. Our 
approach can be extended to other observables in hadronic collisions. 
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I. INTRODUCTION 



The understanding of parton propagation through the nuclear environment is crucial 
for the interpretation of physics phenomena observed at the Relativistic Heavy Ion Col- 
lider (RHIC) and future Large hadron Collider (LHC). The observed strong suppression of 
high transverse momentum hadron at RHIC was considered to be an evidence for the QCD 
quark-gluon plasma The suppression was believed to be the result of medium induced 
radiative energy loss of high energy partons [2|. However, recent data indicate that heavy 
quarks would have to lose the same amount of energy as that of a light quark if the radiative 
energy loss is the only source of the suppression [3[]. On the other hand, we expect heavy 
quarks to lose much less energy than a light quark because of its mass 0]. This discrepancy 
attracted significant theoretical interests in searching for other causes of parton energy loss 
Several studies suggested that the collisional energy loss could be an important source 
of the observed discrepancy [6j]. Because of QCD confinement, we can not observe partons 
directly in experiments, instead, we can only observe final state hadrons. The production 
rate could be affected by the interference of two scattering amplitudes with the same initial 
and final hadronic states but with different partonic interactions. In this paper, we study 
the effect of such quantum interference between two amplitudes with the same as well as 
different multiple parton-level scattering, and show that the interference leads to a suppres- 
sion in single hadron production rate, which might be interpreted as the collisional energy 
loss @- 

The hadronization of partons is a non-perturbative process. However, when the energy 
scale of the scattering Q is much larger than a typical hadronic scale l/fm~ ^ 2 qcdi ^ * s 
the QCD factorization that allows us to separate the calculable short-distance partonic dy- 
namics from the non-perturbative long-distance physics. The effect of the non-perturbative 
hadronization process for a parton of flavor / fragmenting to a hadron h is expressed in 
terms of an universal fragmentation function Df^h(z, Hf), where z is momentum fraction 
of the parton carried by the hadron and /i^ is the fragmentation scale. QCD perturbation 
theory predicts the evolution and the scale dependence on /ip of the fragmentation function. 

Unlike in the vacuum, the fragmenting parton in nuclear medium can have rescattering 
before the formation of final-state hadrons. The effect of the interaction between the nuclear 
medium and the propagating parton will manifest itself as changes in parton fragmentation 
functions. The rescattering can induce extra radiations, and consequently alter the evolution 
of the fragmentation functions 0, [1] , and results in effective parton energy loss 0] . 

However, quantum mechanically, the parton rescattering in a nuclear medium does not 
have to induce radiation. Such rescattering can also alter the production rate of the frag- 
menting parton if it is quantum coherent with the first scattering. Although coherent rescat- 
tering is formally suppressed by additional powers of the hard scale, it could be important if 
the life time of the fragmenting parton is long enough. We show in this paper that quantum 
interference of two scattering amplitudes with different parton-level rescattering reduces the 
production rate of the fragmenting parton, which leads to the suppression of hadron produc- 
tions. In the following sections, we show that such suppression could be effectively expressed 
as a shift in z for the parton fragmentation functions. This result is complementary to the 
radiative parton energy loss induced by the rescattering of the fragmenting parton in the 
nuclear medium. 

In general, the production of hadrons can come from both quarks and gluons in hadronic 
production. However, in the semi- inclusive deep- inelastic lepton-nucleus scattering (SIDIS), 
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FIG. 1: Factorization for semi-inclusive DIS cross section: (a): the leptonic tensor L^ u ; (b): the 
semi- inclusive hadronic tensor. 

the leading hadron production is dominated by quarks. It is the ideal place to study the 
quark energy loss and the knowledge obtained will help us learn more about the gluon energy 
loss in hadronic productions. 

Our paper is organized as follows. In the next section, we first present our derivation of 
double scattering effect. Then in Sec. Ill, we show how to generalize it to the nth-scattering 
and sum up all possible numbers of scatterings, and present the result after the summation. 
In Sec. IV, we extend our results for light quarks to quarks with finite mass and discuss 
the effect of quark mass. In Sec. V, we compare our results with experimental data from 
HERMES. We also discuss the applicable ranges of z value for our result and propose a 
model for larger z region. Finally, in Sec. VI, we summarize our work and discuss extensions 
of our approach to other observables in hadronic collisions. 

II. DOUBLE SCATTERING CONTRIBUTION 

We consider the semi-inclusive DIS production of a single hadron of momentum ph, 



where ki and &2 are the four momenta of the incoming and the outgoing leptons respectively, 
p is the average momentum per nucleon for the nucleus with the atomic number A. To study 
the nuclear effects, we compute the hadron production rate per DIS event, 



where da e A~^ e x / 'dx sdQ 2 is the inclusive DIS cross section. In Eq. fl2]), the Bjorken variable 
xb = Q 2 /{2p ' q) with the virtual photon momentum = (hi — k^ and Q 2 = —q 2 . We 
work in photon-nucleus frame, and choose the target momentum p along F-axis, such that 
pP, _ (jp ,p x ,py ,p z ) = (P, 0,0,P), and only p + = (p° — p z )/\/2 is nonvanishing, neglecting 
target mass. In this frame, the nucleus is moving in the "+" direction. The struck quark 
propagates along the "— " direction and could interacts coherently with the "remnants" of 
the nucleus. The hadron momentum fraction Zh is defined as 
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FIG. 2: Lowest order single scattering Feynman diagram for the semi- inclusive hadronic tensor 
with (a) and without (b) the nonperturbative parton-to-hadron fragmentation attached. 



With the approximation of one-photon exchange, the semi-inclusive DIS cross section 

dc ™ x - 1 e ' zrfaM^, (4) 



dxBdQ 2 dzh 87r x 2 B s 2 Q 2 ' dzh 

where s = (p + ki) 2 is the total invariant mass of the lepton-nucleon system. In Eq. (J3]), the 
leptonic tensor L^ v is given by the diagram in Fig. [Tk, 

L^(kt, k 2 ) = l - Tr( 7 • fc l7 " 7 ■ k 2l v ) . (5) 

The semi-inclusive hadronic tensor for single scattering contribution is represented 
in Fig. [T)d and can be factorized as 

dW^ ^ f 2 f dx 2xp-p h 2xp-p h 
— = £ J dzdz x V^ h {z,z u Q ) J -5(z--^)S( Zl _) 

x<j )q (x,Q 2 )H^(x,z,z 1 ), (6) 

where J2 q runs over (anti) quark flavors, (fi q is the leading twist nuclear quark distribution, 
and Hp, u is the partonic part. In Eq. ([6]) we introduced a two-momentum quark-to-hadron 
fragmentation density V q ^h(z, zi,Q 2 ), which is defined as f dz\b{z — zi)V q _,h(z, %i,Q 2 ) = 
D q ^h{ z i Q 2 ) with D q _,h the normal quark-to-hadron fragmentation function. The in- 
dependence in Eq. (^j is implicit in the argument of the 5-function, 2xp ■ Ph/Q 2 = (x/xB)zh- 
In Fig. [2] we show the lowest order Feynman diagram for the semi-inclusive hadronic tensor 
on a quark state with (a) and without (b) the nonperturbative parton-to-hadron fragmen- 



tation attached. Fig. [2b gives the lowest order Hffi 

H$ = ~e T uu y xe 2 6 (x- -^-) (7) 

where e q is the quark fractional charge, and eF is defined as 

eL = -— \PnQp + QnPu] + ^-PvPv ~ 9/j.v ■ (8) 
^ p ■ q p- q 

Thus, at the lowest order, 

dw!B i 
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< <P q (xB,Q 2 )D q (z h ,Q 2 ), (9) 




FIG. 3: Sample diagram shows the interference between two amplitudes with the same initial and 
final states but different partonic scattering. 



with Zh defined in Eq. (J3j). 

Single scattering is localized and its medium size dependence is limited to that in the 
nuclear parton distribution function. We are interested in additional nuclear effects from 
multiple scattering of the scattered quark and the interference of amplitudes with different 
partonic scatterings, as sketched in Fig. [3j In this section, we first consider the double 
scattering contribution to the semi-inclusive hadronic tensor. 

Fig. H] shows leading order double scattering Feynman diagrams for the semi-inclusive 
hadronic tensor on a quark state. The fermion lines with a short bar represent the contact 
terms of quark propagators, and diagrams with both gluons attached to the incoming quark 
line vanish [12]. Although all diagrams in Fig. H] could contribute to the production rate 
of the hadron in SIDIS, as we discuss below, only diagrams in Figs. Ufa) and (b) can give 
the leading power contribution in the A 1//3 -type medium size enhancement [l3[ . In terms of 
Feynman diagram, contribution of ever y p ropagator consists of two parts: a potential pole 
contribution and a contact contribution 12]. For example, a quark propagator of momentum 
k can always be written as 

ry • k «7 • k ry • n k 2 



k 2 + ie k 2 + ie 2k ■ n k 2 + ie 

where k 2 = and n M is any auxiliary vector with fc-n^O. The first term in the right-hand- 
side of Eq. (11 01) corresponds the potential pole contribution when k 2 — > 0, while the second 
term is the contact contribution [12| . Attaching one gluon to the initial quark line introduces 
a quark propagator, and this propagator will have both the pole and contact contributions. 
The pole contribution is long-distance in nature, representing the interactions between the 
quark and the gluon long before the hard collision between the quark and the virtual photon. 
The pole contribution of the incoming quark propagator should be part of the nuclear quark 
distribution, and is partially responsible for the relatively weak A-dependence of the leading- 



twist parton distributions in a nucleus [14], [15] . On the other hand, the contribution of the 



contact term is localized in space [12j, and does not result into the A 1 / 3 type of nuclear 
enhancement lq . Since Feynman diagrams in Fig. HJ[c) form a gauge invariant subset of 
short-distance partonic contributions, this set of diagrams will not contribute to the leading 
A 1 / 3 -type of nuclear enhancement. Diagrams with the initial-state contact interactions in 
Figs. IU[a) and (b) vanish, while the other two diagrams with final-state rescattering will 
have at least one pole contribution. 

The pole contribution from diagrams with the final-state rescattering, shown in Figs. S(a) 
and (b), is responsible for the leading A 1 ^ 3 -type of nuclear enhancement, because taking the 
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FIG. 4: Leading order double scattering Feynman diagrams that contribute to the semi-inclusive 
deep inelastic scattering on a quark state. Diagrams of (a) and (b) give explicit ^4 1,/3 -type medium 
size enhancement; and diagrams of (c) give localized contributions. 



residue of the unpinched pole effectively sets the gluon momentum to zero and leaves the 



corresponding coordinate space integration of the gluon field to the size of nucleus [13 



Similar to Eq. ([6]), the factorized form of the leading pole contribution from the double 
scattering diagrams in Figs. H](a) and (b) can be expressed as 



= J dz dziV q ^ h (z,z 1 ,Q ) J dxdx 1 dx 1 b{z — 2 — )o\Z\ ^ — , 



x M A (x, xi, an) Hffl(x, fi, an, x B , z, z x ) . (11) 
In Eq. (II ip . the hadronic matrix element M A is given by 



I ! f dy dy 1 djji ^ xp + y - i(x 1 -x) P +y 1 - -i(xi-xi)p+ y ^ 



M A {x,x 1 ,x 1 ) = / — - — e 

X\ — X\ — It X\ — X — le J Z7T Z7T Z7T 

x (P A \m^F +a (y;)F a + (yTMy-)\PA) ■ (12) 

Here we used F +a (v~) = n p d p A a {y~) in light-cone gauge. We adapt the ie prescription 
introduced in Ref. [18| for the two poles l/(an — x\) and l/(an — x). If the partonic part, 
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FIG. 5: Lowest order double scattering Feynman diagrams that contribute to the hard part of the 
leading power ^4 1//3 -type nuclear enhancement. 



Hj£\x, x~i, Xi, Xb, z, Zi), is a nonvanishing smooth function at the poles, which will be veri- 
fied later, the two unpinched poles in Eq. (j!2p can be used to perform the contour integration 
for dx\dx\ in Eq. ffTTj) . In this sense, the hadronic matrix element Ma can be effectively 
written as 

Ma{x, x~i, xi) — > 8(xi — x) 8(xi — fi) J 7 A (x, fij Xi) (13) 
with the function Ta defined as 

F A (x, fi, xi) = J ^ ^ e ixp+y ~ jfr-*)p + v~i- e -i&-*i>p + vl 

x(27r) 2 ^(yr)^(yr)(PA|^(o)^F +Q (yr)^ + (yr)^(^)l^). (u) 

The partonic part Hffl in Eq. (fTTj) is given by the Feynman diagrams in Fig. [5] with the 
quark lines contracted by 57 • p and gluon lines contracted by |<i a/ g. The transverse tensor 
d a p = —g a /3 + n a np + n Q n^ with two lightlike vectors, n = (n + ,n~,nr) = (0,1, 0j_) and 
n = (1, 0, 0_l). The leading pole contribution of Hffi with cut on the left side is 



tt(d) i t / 1 2\ f2na s \ 1 5{x-x B ) 



and the contribution with cut on the right side is 



tt(D) I T f 1 2 1 O^i-XflJ , . 

In-cut = V I -e, ) ( —5— ) t^xb — - — , (16) 



1 2 \ / 2ira s \ 1 £(xi - x 5 ) 
2%/ V 3 J Q* XB x-x x 

with the quark fractional charge e q and e? given in Eq. (jSJ). As we can see, the individual 
contribution of Eqs. ([15]) and (fT6l) is divergent. However, the sum of the contributions is 
finite. 

Summing up contributions from the both cuts, and using Eqs. (fTTp . f[T3"j) . ffTBT) . and ffTBl) . 



we have the leading order contribution from the double final state scattering: 



dW, 
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(17) 



J dx dxi !Fa{x, x, xi) J dz dz\8{x — z-^- 



Ph 
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Q 1 
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x5(z - z x )V q ^ h { z ^ z ^Q 2 ) 



x B Q 2 



Q 2 2p ■ Ph 



S(z - z h ) 5(z 1 - z h ) 



Z\ — z 



Z — Zi 



(18) 



with an overall constant factor (e^/2) e 2 (27ra s /3) from the hard parts in Eqs. (TO)]) and HTB]) 
and a sum over all quark and antiquark flavor. Because of the 8(z — zi), we can expand the 
Z\ in above expression inside the square bracket around z: 



S(z- z h ) 6{zi 



Zh 



Z\ — z 



S'(z-z h ) . 



(19) 



Combining Eqs. ( |T8l) and ( |T9l) . and carrying out all integrations by using the S- functions, 
we have 



dW L 
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T q A (x B ,Q 2 )D q ^ h (z h ,Q 2 



Zh 



q5 T qg (x B , Q 



dD q ^ h (z h , 
dz h 



(20) 



where Yl q runs over all quark and antiquark flavors, and the twist-4 quark-gluon correlation 
function is defined as hj 



T q A g (x B ,Q 2 ) 



dy 

2tt 



ix B p^y 



x(PAi^(o)^-F +a ( y r)^ + (i7r)^(r)l^) . 



dyi dy 1 
2tt 



(21) 



This result shows that interference of amplitudes with different parton-level multiple scat- 
terings affects the hadron production rate in SIDIS. The effect is sensitive to the slope of 
incoming parton flux as well as the shape of the fragmentation function. When combined 
with the lowest order in Eq. Q, the first term in Eq. (1201) is responsible for the high twist 
shadowing of inclusive deep inelastic scattering (DIS) [lOj]. The same effect should also ap- 
pear in the denominator of Eq. ([2]), and therefore, it will not be included in the rest of the 
discussion for the hadron production rate defined in Eq. ([2]). 

From Eq. (fTTj) . we can also derive Eq. (120]) by first expanding x± around x utilizing the 
5-function 5(xi — x): 

5(x — xb) S(xi — xb) 



-5'(x — xb) 



(22) 



X\ — X X — X\ 

and then integrating over dxdx\ and dzdz\. In next section, when we calculate the higher 
order multiple scattering effect, we will follow this derivation to make the presentation 
simpler. 



III. GENERALIZE TO HIGHER ORDER MULTIPLE SCATTERING 



To compute the effect of higher order final state multiple scattering, we add pairs of 
gluon interactions to the struck quark and convert the gluon field operators in the hadronic 
matrix element of W^ u to the corresponding field strength. Each pair of gluon interaction 



will contribute a factor of [10 

dy~ dyf e^ Xi ~ £ ^ p+y i e Kxi-xi-i)p + yf 



x B 



2ira s 

w 



-1 
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i - <'/;-■": v ' (23) 
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Xi—XB—1£ 



"L" ("R") means the gluon pair are to the left (right) of the final state cut line. 

To obtain the leading pole contribution for the partonic part with n additional scattering, 
we need to sum over all diagrams with all possible insertions of the n gluon pairs to both 
sides of the final state cut line. Similar to Eq. ([TBI we can replace the poles in Eq. (|23|) by 
corresponding 5- functions, expand all Xi of 5(xi — x) around x, and obtain 



fJ,V 



2 " 



It: a. 



x B 



(24) 



Convoluting H$ with the fragmentation function and following similar derivations for the 
double scattering case, we obtain the leading pole contribution for semi-inclusive hadronic 
tensor, with n additional scattering: 



dW$ 

dz h 



-e T V, 



z h - 



-M^(x B ,Q 



dzt 
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(25) 



with the multi-field matrix element M\ given by 



M^(x,Q 2 



dy, 



2tt 



^ e ^o-(p A |^( )x_v/(%) n 
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i=l 



P + dyre(yr)F\yr) 



\Pa)- (26) 



The integration J p + dy i in Eq. (1261) gives the nuclear size dependence 15|. And the operator 
F 2 (yl) is given by 



dy~ F+^yDF+iyi 



(27) 



2tt p+ 

Compare with the operator definition of gluon density, we can see that its expecta- 
tion value can be related to the small-x limit of the gluon distribution, (p\F 2 (y~)\p) w 
lim^o | x G(x, Q 2 ), and is independent of yi Tol . 

In order to evaluate the multi- field matrix element in Eq. (|26|) . we approximate the 
expectation value of the product of operators to be a product of expectation values of the 
basic operator units in a nucleon state of momentum p = Pa/ 'A: 

n n 

(Pa\ O H 6, \P A ) = A (p | 6 I p) n [ N p (P \Oi\p) 



1=1 



i=l 



where N p is the normalization. In a model of constant lab frame nucleus density p(r) 
3/(47rrg), we have 



P + dyr9(yi)N p (p\F 2 (yr)\p) 



9 



lQnr 2 



(A^-l)(p\F 2 (yn\ P ) 



(28) 



The factor (A 1 / 3 — 1) is taken such that the nuclear effect vanishes for A = 1. With the 



above model for \ we have 
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The quantity n 2 represent the characteristic scale of quark interaction with the medium lol 
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Summing the y4 1,/3 -enhanced contributions in Eq. f[2"9~l) to all order in n, we have 
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Q 2 ' W 



(31) 



where N is the upper limit on the number of quark-nucleon interactions. In deriving Eqs. (I3T|) 
we have taken N pa oo because the effective value of k 2 is relatively small. Eq. (|3T|) is the 
main result of this paper. It shows that the net effect of multiple scattering without induced 
radiation for a propagating quark in the medium is equivalent to a shift in the variable z 
for the quark fragmentation function D q ^h{z), which leads to a suppression of the hadron 
production rate. This is a result of the quantum interference of amplitudes of multiple 
scattering. Such a shift in z for the fragmentation function is very similar to the effect of 
the parton energy loss model proposed in Ref. [21]. The shift 



Az(z h ) = z h 



ft 2 (A 1/3 - r 

Q 2 



(32) 



oc 



depends on only one parameter k 2 and the medium length. The parameter 
lim^o xG(x, Q 2 ) with G(x,Q 2 ) the gluon distribution function. It can be related to the 



A 2 in LQS model for twist-4 quark-gluon correlation function T qg (x) = A A /' s \ 2 (j) q {x) 
by k 2 = (47r 2 a s /3)X 2 . The A 2 has been estimated using Drell-Yan transverse momentum 
broadening and DIS momentum imbalance 17, 2(|, and was in the range of 0.01 —0.1 GeV 2 . 

Using our result given in Eq. (13T1) ,we obtain the hadron production rate defined in Eq. (J2J) 
for SIDIS on a nucleus target A: 



E 



e 2 6 A 



(x B , Q 2 ) D q ^ h (z h + Az(z h )) 



E 



K X B,Q 2 ) 



(33) 



with Az(z h ) given in Eq. (|32|) . 



IV. THE QUARK MASS EFFECT 



In the above derivation, we concentrated on the multiple scattering of light quarks and 
ignored the quark mass. If the stuck quark is a heavy quark, we then can not ignore the 



quark mass. In this case, the initial quark momentum can not be approximated as xp. 
Instead, we have the quark momentum 



W 1 = xp + n^ + 



m 



2xp^ 



-n r 



(34) 



with m the quark mass. Due to the on-shell condition (k+q) 2 = m 2 , the final state 5-function 
is modified as 

S ((k + qf - m 2 ) = = == ^ 5(x - x B ) (35) 



with 
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In addition, for additional scattering with each gluon pair in the medium, the interaction 
will contribute a factor 



XBr, 
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Correspondingly, due to the multiple scattering with the medium, the shift Az(zh) in frag- 
mentation function for a heavy quark with mass m is 



Az m (z h ) 



1 



'l + 4m 2 /Q 2 



k 2 (A 1/3 - 1) 

Q 2 



(38) 



From Eq. ( |38j) , we see that the z-shift for a heavy quark has a similar functional form to that 
of a light quark, except that it has an extra factor that depends on the quark mass. This 
mass dependent factor makes the z-shift smaller for a heavier quark. However, as we can 
see from Figj6] and Figd the fragmentation function of a heavy quark have very different 
z-dependence from that of a light quark. The z-shift in fragmentation function can result 
in very different hadron production ratio, even though they have similar forms of z-shift. 
Therefore, we expect that the net effect of multiple scattering for heavy meson production 
to be very different from that for light mesons. 

V. NUMERICAL ESTIMATES AND COMPARISION WITH DATA 



In order to compare with data [22], we compute the ratio of R A for a nuclear target A to 
that of a deuterium target D: 



R 



M 



R D 

V e 2 <pf(x B , Q 2 ) D q ^ h (z h + Az(z h )) V e 2 tf?(x B , Q 2 ) 



E Q e 2 ^(x B ,Q 2 )D^ h (z) E„e^(* B ,Q 2 ) " 



(39) 




FIG. 6: u-quark fragmentation function given by Ref. [251 ] . 




The super script "A" and "D" represent the nuclear target of atomic weight A and the 
deuterium target, respectively. 

To obtain the numerical estimate of the double ratio Rm in Eq. fl39|) . we use the lowest 

For the nuclear dependence of parton distribution, 
Figs. IBlfTTl compares our result with the 



order CTEQ6 parton distributions [23 



we use the parameterizations given in Ref. [24 



data from HERMES experiment 22j] . The experiment data points have xb in the range of 
0.084 - 0.1, and the Q 2 ~ 2.2 - 2.6 GeV 2 . In plotting Figs. [8lfTTl we used the fragmentation 
function provided by Ref. 



25 



The dot lines represent the result with double scattering 
only. The dashed curves represent the result when we sum to all order. We can see that 
our curve is slightly above the data, because our calculation only include multiple scattering 
without induced radiation. The induced radiation will give further suppression 0, @] and 
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FIG. 8: The double ratio Rm of Eq. (|5U|) for tt + production with Krypton target, compared with 
HERMES data. 
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FIG. 9: The double ratio Rm of Eq. (|39p for tt production with Krypton target, compared with 
HERMES data. 



brings down the curve. We also notice that at large Zh region, our curve is steeper than the 
data point. This is because our result sums the additional scattering to all order and it is not 
applicable for large Zh- At large Zh, the hadron forms early, and have shorter formation time. 
In this case, Summing the additional scatterings to all order is an unrealistic approximation. 
In order to take into account of the hadron formation time, which is proportional to (1 — Zh) 




FIG. 10: The double ratio Rm of Eq. (f5Uj) for K + production with Krypton target, compared with 
HERMES data. 



26j . we modify the Zh shift to be at large Zh 

. , , K 2 (A^ - 1) (1 - z h ) 

Az(z h ) = z h V '- V _ h! when z > z c . (40) 

iqj 1 Zq 

In Eq. ( HO"]) . z c is a parameter. In this model, we assume that when z^, < z c , the hadrons 
will form outside of the nucleus, and we do not need to worry about the hadron formation 
time. We can apply our resumed result of Eq. (|32|) when z h < z c . When z h > z c , the hadron 
formation process may start early, and we use the modified Az(zh) given in Eq. (140!) to take 
into account of the formation time. The solid lines in Figs. IBlfTTl show our estimates when we 
use the above modified Az, and choose z c = 0.6. The curve is above the data points at large 
Zh, because here we did not consider the nuclear absorption of the pre-hadron state. The 
nuclear absorption of pre-hadron state should give additional suppression and bring down 
the curve a little more (26[. 

In Fig. [121 we plot the figure for D meson production in eA scattering. Due to the 
limit in collision energy, there is no data for D meson production from Hermes experiment. 
However, future Electron-ion Collider (EIC) experiments at Brookhaven should be able to 
observe the semi-inclusive D meson production [27|]. In obtaining Fig. [121 we used the 
c-quark fragmentation function of Ref. j28|. From Eq.(39), we see that the ratio Rm is 
mainly determined by the shift Az(zh) and the shape of fragmentation functions. Since 
the fragmentation function for light- and heavy-meson have very different z-dependence, as 
shown in FigEJand FigJTJ the shift in Zh results in very different double ratio Rm- At smaller 
Zh region, Zh may actually result in Rm > 1 for D mesons, due to the characteristic shape 
of the c-quark fragmentation function, as we can see from Fig. [71 

In Fig. l8lfT2l we only illustrate the size of the suppression due to the shift of Zh, which is 
caused by coherent multiple scattering without induced radiation. A more complete analysis 
should include other effects, such as induced radiation. 




FIG. 11: The double ratio Rm of Eq. (f5Uj) for K production with Krypton target, compared with 
HERMES data. 
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FIG. 12: Our estimate of double ratio Rm for D meson production with gold target at Q 2 = 9 
GeV 2 and Q 2 = 25 GeV 2 . 



VI. SUMMARY AND OUTLOOK 

In summary, we demonstrated that the coherent multiple scattering of a propagating 
quark with the medium, without the induced radiation, can also change the quark frag- 
mentation or hadron's production rate. The net effect of leading power contributions in 
medium length is equivalent to a shift in the fragmentation function's Zh — * Zh + A2. At 



the lowest order, the shift Az is given by an universal matrix element, which is proportional 
to lim x ^oxG(x, Q 2 ). We also show that for a quark with mass m, the shift will be smaller. 
Our result could be interpreted as the collisional energy loss, which is complementary to the 
energy loss of induced radiation. However, beyond the leading order, the separation of the 
collisional energy loss and that of induced radiation will depend on the factorization scheme 
and will not be unique, and need further study. 

In this paper, we derived the effect of the radiationless multiple scatterring in semi- 
inclusive DIS, our approach can be systematically generalized to hadron production in p + A 
and A + A collisions. Because of the convolution of two parton distributions, which are 
steep falling functions of parton momentum fraction x, the hadron production in hadronic 
collisions is dominated by the large z part of fragmentation functions, in particular, for 
heavy meson production [29| . As shown in Fig. 0, the heavy quark fragmentation functions 
are steep falling functions of z in large z region. As a result, depending on the momentum 
of the observed heavy meson, or the effective range of z, we expect that the radiationless 
multiple scattering will lead to a suppression of heavy meson production, similar to that in 
light hadron production. 
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